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replication in human airways, the OC43 
hemagglutinin-esterase lost its receptor¬ 
binding function. Consequently, virion- 
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SUMMARY 

Human betal-coronavirus (pICoV) OC43 emerged 
relatively recently through a single zoonotic introduc¬ 
tion. Like related animal pi CoVs, OC43 uses 9-O-acet- 
ylated sialic acid as receptor determinant. PICoV 
receptor binding is typically controlled by attach¬ 
ment/fusion spike protein S and receptor-binding/re- 
ceptor-destroying hemagglutinin-esterase protein 
HE. We show that following OC43’s introduction into 
humans, HE-mediated receptor binding was selected 
against and ultimately lost through progressive accu¬ 
mulation of mutations in the HE lectin domain. Conse¬ 
quently, virion-associated receptor-destroying activ¬ 
ity toward multivalent glycoconjugates was reduced 
and altered such that some clustered receptor popula¬ 
tions are no longer cleaved. Loss of HE lectin function 
was also observed for another respiratory human co- 
ronavirus, HKU1. This thus appears to be an adapta¬ 
tion to the sialoglycome of the human respiratory tract 
and for replication in human airways. The findings sug¬ 
gest that the dynamics of virion-glycan interactions 
contribute to host tropism. Our observations are rele¬ 
vant also to other human respiratory viruses of zoo¬ 
notic origin, particularly influenza A virus. 

INTRODUCTION 

Coronaviruses (CoVs), long considered of high veterinary impact 
exclusively, are now generally recognized as zoonotic threats of 
pandemic potential in consequence of the 2002/2003 outbreak 
of severe acute respiratory syndrome (SARS) and the emer¬ 
gence of Middle East respiratory syndrome (MERS) in 2012 (de 
Wit et al., 201 6). SARS-CoV was contained within 3 months after 
its discovery, while MERS-CoV causes a classical zoonotic 
infection with limited human-to-human spread and, as yet, is 
incapable of sustained community transmission (Reusken 
et al., 2016; Zumla et al., 2015). However, four other respiratory 


coronaviruses—alphacoronaviruses NL63 and 229E and lineage 
A betacoronaviruses OC43 and HKU1 —successfully breached 
the species barrier and are currently maintained in the human 
population worldwide through continuous circulation (Su et al., 
2016; de Groot et al., 2011). Conceivably, the study of these 
genuine human coronaviruses (HCoVs) may yield clues to what 
is required for viral adaptation to the human host and thereby in¬ 
crease our understanding of the probabilities and risks of coro- 
navirus cross-species transmission. 

OC43 and HKU1, while generally associated with benign com¬ 
mon colds in healthy immunocompetent individuals, may cause 
significant morbidity and even mortality in the frail (Morfopoulou 
et al., 2016; Woo et al., 2005a). OC43, the best-studied HCoV, 
apparently arose relatively recently, 120 to 70 years ago, with 
the most recent common ancestor of all extant OC43 variants 
dating to the 1950s (Lau et al., 2011; Vijgen et al., 2005, 2006). 
OC43 groups in the species Betacoronavirus-1 (pICoV), 
together with highly related viruses from ruminants (bovine coro- 
navirus, BCoV), swine (porcine hemagglutinating encephalomy¬ 
elitis virus, PHEV), equines (equine coronavirus, ECoV), leporids 
(rabbit coronavirus HKU14, RbCoV), and canines (canine respi¬ 
ratory coronavirus, CRCoV) (de Groot et al., 2011; : igure 1A). 
The extraordinary radiation of pi CoVs might be explained from 
their receptor usage, as they attach to 9-O-acetylated sialic 
acids (9-O-Ac-Sias) (Matrosovich et al., 2015), i.e., glycan com¬ 
ponents common in mammals and birds (Traving and Schauer, 
1998). Paradoxically, however, most pi CoVs, including OC43, 
have very narrow host ranges. They form distinct monophyletic 
clades congruent with host selectivity, and phylogenetic evi¬ 
dence strongly argues against recurrent inter-species transmis¬ 
sions (Lau et al., 2011; Vijgen et al., 2006; see also Figure 1A). 
With the emergence of OC43 seemingly sparked by a singular 
one-time zoonotic event, the founder virus likely possessed 
unique traits already to allow efficient infection of, and transmis¬ 
sion among, humans. In turn, these traits and subsequent 
adaptations to the new niche must have closed the door on rein¬ 
troduction of the virus into animals. Despite the close genetic 
relationship between OC43 and BCoV-Mebus, particularly 
(96.6% identity across their entire genomes), their high-preva¬ 
lence endemicity/enzooticity, and the scale and frequency of 
interactions between their host species, there is no evidence of 

CrossMark 


356 Cell Host & Microbe 21 , 356-366, March 8, 2017 © 2017 Elsevier Inc. 





CelPress 


PICoV 



B 



• 

BCoV 

• 

OC43 USA/1967 

□ 

CRCoV 

V 

OC43 NL/A/2005 

♦ 

PHEV 

A 

OC43 NL/B/2005 

▲ 

ECoV 

♦ 

HCoV-HKUl 

♦ 

RbCoV 




C 


Standard HAA 


BCoV 
CRCoV 
RbCoV 
OC43 USA/1967 
OC43 NL/A/2005 
OC43 NL/B/2005 
PHEV 
ECoV 
HKU1 


oooooooo-e- 

OOOOOO O. v 
OOOOOOOCK 


O >00000®'®®®©‘ 

OOOOOOOOO 6® 

. . . 



\ r 



HE°-Fc 


q Nanobead HAA 

BCoV oooooooooooo 

PHEV OOOOOOOOOOOO 

BCoV F 211 A GO'.' • V s V 

OC43 USA/1967 ooooo 
OC43 NL/A/2005 fcfO 
OC43 NL/B/2005 
HKU1 

Beads only CO '.' 



HE°-Fc 


Figure 1. Loss of HE-Mediated Receptor 
Binding in Human Betacoronaviruses 

(A) Evolutionary relationships among lineage A be¬ 
tacoronaviruses. Neighbor-joining tree based on 
(3CoV lineage A ORFIb sequences in the NCBI 
database (n = 206), with 100% bootstrap support 
for all major branches. Evolutionary distances were 
computed using the Maximum Composite Likeli¬ 
hood method in MEGA6 (Tamura et al., 2013). The 
positions of human coronaviruses OC43 and 
HKU1 are highlighted relative to those of animal 
lineage A betacoronaviruses, including the various 
(31CoV subspecies and classical mouse hepatitis 
virus-type murine coronavirus (MuCoV). ChRCoV- 
HKU24, Chinese rat coronavirus HKU24 (Lau et al., 
2015); LAMV, Longquan Aa mouse coronavirus 
(Wang et al., 2015); LRLV, Longquan Rl rat coro¬ 
navirus (Wang et al., 201 5). The position of CRCoV, 
a recent split-off of BCoV, is indicated by a red dot. 
(See also r igure SI .) 

(B) HE°-Fc lectins (2-fold serial dilutions, starting at 
50 ng/|.iL) were compared by sp-LBA for relative 
binding to BSM (with 50 ng/fiL BCoV-HE°-Fc set 
at 100%). 

(C) Conventional HAA with 2-fold serial dilutions of 
HE°-Fc fusion proteins of (31CoV members and 
HCoV-HKUl (starting at 25 ng/well). Wells positive 
for hemagglutination are encircled. 

(D) High-sensitivity nanobead HAA. Non-com- 
plexed nanobeads (“beads only”) and nanobeads 
complexed with lectin-inactive mutant BCoV 
HE 0 F 211 A were included as negative controls. 


OC43 spreading to cattle nor of BCoV (or any other pICoV) 
efficiently spreading into humans ( : igure 1 A). 

pICoV attachment involves two distinct types of surface pro¬ 
jections: large 20-nm “spikes” comprised of homotrimers of the 
class I fusion protein S and stubby 8-nm protrusions that are ho- 
modimeric assemblies of the hemagglutinin-esterase protein HE 
(de Groot et al., 2011). The S protein, common to all corona- 
viruses, mediates receptor binding and fusion of the viral and 
host membranes (Heald-Sargent and Gallagher, 2012). HEs are 
found in toroviruses and in influenza C and D viruses, but among 
coronaviruses, only in lineage A betacoronaviruses (de Groot, 
2006; de Groot et al., 2011; Hause et al., 2014; Matrosovich 
et al., 2015). HE monomers have a bimodular structure with a 
carbohydrate-binding (“lectin”) domain appended to an enzy¬ 
matically active sialate-O-acetylesterase (“esterase”) domain 
(Langereis et al., 2009; Rosenthal et al., 1998; Zeng et al., 
2008). Typically, in pICoVs, both the spikes and HEs bind 
9-O-Ac-Sias, while the HE esterase domain promotes virus 
elution through receptor destruction. In consequence of the 
opposing activities of receptor binding and receptor destruction, 
pICoV attachment to sialoglycans is dynamic and reversible. 
Thus, dead-end binding of virions to decoy receptors in oropha¬ 
ryngeal, respiratory, and gastrointestinal mucus may be pre¬ 
vented. Moreover, in infected host cells, HE-mediated receptor 
destruction is essential for efficient release of viral progeny (Des- 
forges et al., 2013). 

Here we present a comprehensive structure-function study of 
OC43 HE. We demonstrate that over decades after OC43’s intro¬ 


duction, its evolution was marked by a progressive loss of HE re¬ 
ceptor-binding activity through the accumulation of select muta¬ 
tions in the HE lectin domain. The effect of these mutations on 
the organization of the carbohydrate-binding site (CBS) and on 
receptor binding is explained from the crystal structure of the 
BCoV HE-receptor complex (Zeng et al., 2008) and visualized 
directly by the structure of a contemporary OC43 HE solved to 

o 

2.45A resolution. Evidence is provided that loss of HE recep¬ 
tor-binding activity resulted in a reduction of virion-associated 
esterase activity toward multivalent clustered substrates. We 
propose that inactivation of the HE lectin domain altered the 
balance between virion binding and esterase-mediated virion 
elution, apparently as an adaptation to the sialoglycome of the 
human respiratory tract, and we speculate that this may be a 
contributing factor to host selectivity. This view is supported by 
our observation that the HE of HCoV-HKUl also lost its recep¬ 
tor-binding properties. A mechanism is proposed in which 
accessibility of receptors to destruction depends on HE lectin 
function in relation to S-HE size differences and in which the bal¬ 
ance between attachment and catalysis-driven virion elution is a 
determinant of host tropism. 

RESULTS AND DISCUSSION 

Lectin Properties of (31CoV HEs 

A comprehensive set of pICoV HEs, expressed as esterase- 
inactive Fc-fusion proteins (HE°-Fc), was compared for their 
Sia-binding properties. In accordance with previous findings 
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(Langereis et al., 201 5), the HEs of most animal (31 CoVs bound to 
bovine submaxillary mucin (BSM) in a 9-O-Ac-Sia-dependent 
fashion in solid-phase lectin-binding assays (sp-LBA; : igure 1 B). 
Remarkably, however, those of PHEV strain VW572 and proto¬ 
type OC43 strain USA/1967 (also known as ATCC-VR759; McIn¬ 
tosh et al., 196 ) showed no detectable binding. In hemaggluti¬ 
nation assays (HAA), more sensitive than sp-LBA, PHEV HE 0 
tested positive, albeit weakly ( igure 1C). OC43 USA/1967 
HE 0 , however, did not hemagglutinate. To augment HAA sensi¬ 
tivity through multivalency-driven high-avidity binding, we com- 
plexed HE°-Fc chimeras to protein A-coated nanobeads ( ig¬ 
ure ID). For PHEV HE 0 , sensitivity was increased 250-fold 
as compared to the standard assay. Under these condi¬ 
tions, modest, but reproducible, hemagglutination by OC43 
USA/1967 HE 0 was detected ( : igure ID). Apparently, it has 
lost most, though not all, of its lectin function. 

Our observations prompted the question of whether loss of HE 
lectin activity is a strain-specific trait resulting from adaptation to 
in vitro propagation or a characteristic also shared by naturally 
occurring OC43 viruses. We therefore RT-PCR-amplified HE 
genes from more recent sputum-derived OC43 strains (respira¬ 
tory season 2005). The encoded proteins OC43/NL/A/2005 
and OC43/NL/B/2005 HE 0 , representative for the two major HE 
lineages in OC43 (designated “A” and “B,” Figure SI and vide 
infra), did not show any detectable binding by sp-LBA or 
nanobead HAA ( igures 1 B-1D). Apparently, loss of affinity for 
9-O-Ac-Sias as seen for OC43/USA/1967 HE is not an artifact. 
In fact, the data suggest that HEs of contemporary OC43 vari¬ 
ants have lost receptor-binding activity altogether. 

Lectin Function of OC43 HE Impeded by a Combination 
of Mutations 

The HE ectodomains of BCoV strain Mebus and OC43 
USA/1967, each 365 residues in length, differ at 18 positions 
only figures 2A and 2B). To identify the differences accounting 
for loss of lectin activity in OC43 USA/1967 HE, we introduced 
OC43-specific substitutions in BCoV-Mebus HE°-Fc, in sets 
and individually, and tested for a reduction in lectin activity by 
HAA ( : igure S2A). In a complementary approach, OC43 HE res¬ 
idues were systematically replaced by BCoV HE orthologs ( ig¬ 
ures S2B-S2D). The results show that the mutations in the 
esterase and MP domains of OC43 USA/1967 HE ( igures 2A- 
2C; : igure S2A) do not affect ligand binding. In fact, loss of 
9-O-Ac-Sia binding can be attributed to four out of eight substi- 
tutions in the lectin domain (T 114 N, R 177 P, E 178 Q, and F 247 L; Fig- 
ures S2A-S2D). Combined replacement of these residues in 
OC43 USA/1967 HE 0 by BCoV orthologs restored binding affinity 
to that of BCoV-Mebus HE ( : igure 2C). Individual replacement of 
these residues in the context of a restored OC43 USA/1967 HE 
(“TREF”) showed that each mutation affects receptor binding 
to more (T 114 N, E 178 Q) or lesser (R 177 P, F 247 L) extent (Figure 2C). 

The impact of these mutations on 9-O-Ac-Sia binding can be 
understood from the crystal structure of the BCoV-Mebus 
HE-receptor complex ( !eng et al., 2008). The BCoV HE recep¬ 
tor-binding region is comprised of six surface loops, five grafted 
on the beta-sandwich core of the lectin domain and one 
emanating from the esterase domain ( igure S3A). The actual 
carbohydrate-binding site consists of a deep hydrophobic 
pocket PI and a more shallow hydrophobic depression P2 


that accommodate the methyl groups of the Sia-9-O- and 
-5-A/-acetyl moieties, respectively ( igure 4E). The side chain 
of F 211 (a residue in the pi 2/(313 (3-hairpin) separates PI from 
P2 and, in the HE-receptor complex, intercalates between the 
Sia-acetyl groups. Ligand binding, thus largely based on shape 
complementarity and hydrophobic interactions, is stabilized by 
extensive protein-sugar hydrogen bonding, particularly involving 
the (37—pi 0 p-hairpin, through main chain atoms of L 212 and 
N 214 to the sialate-5-A/-acyl amide and the Sia carboxylate, 
respectively, and through the side-chain hydroxyl of S 213 to the 
sialate-C8-hydroxyl and -C9-acyl oxygen. Y 184 in the p5-p6 
loop (residues 176-185) is particularly important for ligand bind¬ 
ing as its side-chain hydroxyl group hydrogen bonds with the 
sialate-9-O-acetyl carbonyl, while its aromatic ring together 
with the side chains of F 211 , L 266 , and L 267 walls the PI pocket 
that is key to receptor recognition ( -igure 4E). 

Residues of the a5-a6 loop (residues 111 -118) do not directly 
contact the ligand. However, the T 114 N substitution in OC43 
HE created a novel A/-linked glycosylation site (TTS^ NRS) fig¬ 
ure 2A). A glycan attached to this site would extend into the CBS 
and potentially reduce receptor binding through steric hin¬ 
drance. Indeed, the introduction of this glycosylation site in 
BCoV HE 0 resulted in an almost complete loss of binding as 
detected by sp-LBA ( : igure 2C) and a 250-fold loss of apparent 
binding affinity as measured by HAA ( 'igure 2D; = igure S2A). 
The same mutant protein, but now expressed in A/-acetyl 
glucosamine transferase I (GnTI)-deficient HEK293S GnTI - cells 
instead of HEK293T cells, showed a less pronounced loss of 
binding affinity (30-fold in HAA; igure 2D). The data show that 
glycosylation at N 114 indeed reduces the affinity of the HE lectin 
domain for 9-O-Ac-Sias and that this reduction is positively 
correlated with glycan size/complexity. 

In BCoV HE, E 178 fixates the p5-p6 loop by engaging in double 
hydrogen bonds with the backbone nitrogen atoms of S 155 and 
A 156 in the p4-p5 loop ( z igure S3B). We offer that in OC43 HE, 
the replacement by Gin and consequential loss of inter-loop 
hydrogen bonding destabilizes the p5-p6 loop and that its 
increased flexibility impairs the Sia-binding site, presumably by 
affecting the critical PI pocket. The consequences of the other 
two substitutions in OC43 USA/1967 HE are less evident, 
although it can be envisaged that the R 177 P replacement might 
again affect folding and/or stability of the (35-J36 loop and 
thereby indirectly the positioning of the critical Y 184 residue. 
The F 247 L substitution is more difficult to explain as the affected 
residue locates to the core of the HE lectin domain, distal to the 
CBS, and the mutation apparently decreases ligand binding 
through indirect long-range effects. The combined data show 
for the early OC43 isolate USA/1967 that the reduced lectin ac¬ 
tivity of its HE resulted not from a single but from a combination 
of mutations. 

Progressive Loss of HE Lectin Function during OC43 
Evolution 

In the course of this study, many full genome sequences from 
OC43 field variants from the US, Western Europe, and China, 
well documented with respect to place and date of virus sam¬ 
pling, became available in GenBank. This wealth of information 
allowed us to study the recent evolution of the OC43 HE protein. 
In phylograms, in accordance with the presumptive date of 
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Figure 2. Loss of OC43 USA/1967 HE 0 Lectin Affinity Attributed to a Combination of Four Mutations 

(A) Alignment of BCoV, OC43 USA/1967, OC43 NL/A/2005, and OC43 NL/B/2005 HE with sequences color coded by domain (membrane-proximal domain, red; 
esterase domain, green; lectin domain, blue). Residues crucial for esterase activity (SGNDH) are annotated. Amino acid differences are marked in black. 
N 114 substituting for Thr is colored in cyan, and the resulting A/-linked glycosylation site (NRS) marked with red dots. D 220 is marked with a red asterisk. 

(B) Overall structure of BCoV HE (PDB: 3CL5) with mutations in OC43 USA/1967 HE indicated. Domain coloring as in (A). 

(C) Comparison of binding affinities of BCoV HE 0 , OC43 USA/1967 HE 0 , and derivatives by sp-LBA. BCoV, BCoV-Mebus HE 0 ; OC43TREF, OC43 USA/1967 HE 0 
with N 114 , P 177 , G 178 and L 247 replaced by BCoV-Mebus HE orthologs; T 114 N, R 177 P, E 178 G, and F 247 L, TREF derivatives with indicated residues re-converted to 
the autologous OC43 orthologs. 

(D) Conventional HAA with BCoV-Mebus HE°T 114 N expressed in HEK293T or HEK293S GnTr cells. (See also : igure S2.) 
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Figure 3. Complete Loss of HE Lectin Func¬ 
tion during OC43 Evolution due to Progres¬ 
sive Accumulation of Mutations 

(A) Close up of the BCoV-Mebus HE CBS in 
complex with a-Neu5,9Ac 2 2Me (in sticks). Resi¬ 
dues comprising the metal-binding site (MBS) are 
also shown in sticks, the potassium ion is shown 
as a magenta sphere, and interactions between 
the metal ion and coordinating amino acid resi¬ 
dues as red dashed lines. 

(B) Disruption of the MBS in BCoV HE leads to loss 
of lectin affinity. sp-LBA as in Figure 1 B. 

(C) Disruption of the BCoV HE lectin MBS 
through D 220 Y substitution renders receptor 
binding non-detectable even by high-sensitivity 
nanobead HAA. 

(D) The Y 220 D mutation partially restores lectin 
affinity of OC43 NL/B/2005 HE 0 . 

(E) D 220 H substitution in BCoV HE results in ther¬ 
molability of receptor binding. Conventional 
HAA before (4°C) and after (37°C) a temperature 
shift up. 

(F) Receptor binding of OC43 NL/A/2005 HE 0 and 
derivatives as determined by sp-LBA. In OC43 
NL/A/2005 HE 0 , N 114 , P 177 , Q 178 , and L 247 were 
replaced by BCoV-Mebus HE orthologs (TREF), in 
combination with (1) H 220 D substitution (TREF, 
H 220 D), (2) repair of the (35-(36 loop by substituting 

D 183 YYY 186 for ||J (TREFj _|| T)j or (3) |_|220 D and 

repair of the 35-36 loop (TREF, H 220 D, -NT). 
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OC43 emergence, the HE of the USA/1967 strain is placed close 
to the root ( : igure SI). In the decades thereafter, the HE proteins 
divided into three clades. Clade C, the least populous, is 
comprised of HEs from OC43 strains sampled in the early 
1990s with no recent representatives and may have been re¬ 
placed by clades A and B. Comparative sequence analysis in 
combination with structure-function analysis showed that both 
type A and B HEs acquired several mutations in the lectin domain 
additional to those already present in OC43 USA/1967, which 
explains the loss of HE receptor binding in extant OC43 variants. 
Unexpectedly, in B-type HEs, but not in those of types A and C, 
glycosylation at N 114 was lost again due to a S 116 A substitution 
(NRS->NRA). Further inspection revealed the replacement in 
B-type OC43 HEs of D 220 by Y. D 220 is part of a potassium ion¬ 
coordinating metal-binding site (MBS), a signature element of 
coronavirus HE lectin domains involved in the organization of 
the Sia-binding site ( igure 3A) ( [eng et al., 2008). The MBS is 
formed by main-chain carbonyl atoms of S 221 , E 265 , and L 267 , 
together with the side chains of Q 222 , S 263 , and D 220 , the latter 
of which is particularly important as it balances the positive 


charge of the metal ion. In BCoV HE 0 , 
the D 220 Y mutation resulted in loss of all 
detectable binding to 9-O-Ac-Sia ( = ig- 
ures 3B and 3C). In turn, a Y 220 D back 
mutation in the type B HE of OC43 
NL/B/2005 restored Sia-binding activity 
to levels detectable by high-sensitivity 
nanobead HAA ( igure 3D). Again, the 
structure of the BCoV HE-receptor com¬ 
plex provides an explanation for the 
detrimental effect of the mutation. As L 266 and L 267 are in a 
loop that is stabilized through the coordination of the potassium 
ion, disruption of the MBS might affect the folding of this loop, 
the positioning of the Leu residues, and thereby the structure 
of receptor-binding pocket PI. 

Remarkably, also in A-type OC43 HEs, the MBS was targeted 
by substitution of D 220 , now by His ( igure 2A). As evident from 
comparative sequence and phylogenetic analysis, the substitu¬ 
tions in A- and B-type HEs were independent events, and their 
selection must be ascribed to convergent evolution. Again, 
D 220 H mutation in BCoV HE resulted in loss of detectable Sia 
binding as measured by sp-LBA ( igure 3B). Unexpectedly, 
however, in conventional HAA routinely performed at 4°C, 
wild-type and mutant BCoV HEs hemagglutinated to similar 
extent ( : igure 3E). The effect of the mutation became evident af¬ 
ter a shift up to 37°C, upon which hemagglutination with wild- 
type BCoV HE remained stable while that of the mutant protein 
resolved ( : igure 3E). Apparently, substitution of D 220 by His is 
less disruptive than by Tyr, but results in thermolability of HE re¬ 
ceptor binding, likely to cause reduced binding affinity under 
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Figure 4. Crystal Structure of OC43 
NL/A/2005 HE 

(A and B) Side-by-side cartoon representations of 
the overall crystal structures of BCoV-Mebus HE 
(A) and OC43 NL/A/2005 HE (B). HE monomers are 
colored gray or by domain (as in igure 2A). 

(C) BCoV and OC43 HE have identical esterase 
catalytic sites. Overlay of BCoV (blue) and OC43 
NL/A/2005 (green) HE esterase domains. Cartoon 
representations with residues crucial for activity 
indicated as sticks. 

(D) Overlay of BCoV and OC43 NL/A/2005 HE 
lectin domains. F 211 indicated with sticks. 

(E) 9-0-Ac-Sia (with carbon atoms in cyan) binding 
in the BCoV-Mebus HE lectin CBS as observed in 
the crystal complex (PDB: 3CL5). Close up with 
contacting amino acid side chains shown in stick 
representation. Hydrogen bonds are shown as 
black dashed lines, and hydrophobic interactions 
with the Sia-9-0- and -5-/V-methyl groups as thin 
gray lines. PI and P2 indicate the pocket and hy¬ 
drophobic depression, which accommodate the 
methyl groups of the Sia-9-0- and -5-A/-acetyls, 
respectively. 

(F) Close up of the inactivated CBS of OC43 
NL/A/2005 HE. Residues corresponding to those 
in (E) are in stick representation and colored by 
atom type. NAG, A/-acetylglucosamine attached to 
N 114 , is shown in stick representation. (See also 
Figure S3 and 'able SI .) 



physiological conditions. Repair of the OC43-specific substitu¬ 
tions already present in the HE of the 1967 strain ( r igure 2A), in 
conjunction with a H 220 D back substitution, did not restore lectin 
activity of OC43 NL/A/2005 HE ( : igure 3F), suggesting the pres¬ 
ence of at least one additional mutation to prevent receptor bind¬ 
ing. Comparative sequence analysis revealed a major change 
uniformly shared by A-type HEs again in the p5-p6 loop, with 
four adjacent residues 183-186, Asp-Tyr-Tyr-Tyr, replaced by 
lle-lle-Thr ( : igure 2A), apparently as a result of a double frame- 
shift mutation ( igure S3C). Repair of this mutation, in combina¬ 
tion with restorative changes at the other five sites, was required 
for OC43 NL/A/2005 HE to regain receptor-binding activity ( ig¬ 
ure 3F). The impact of the p5-p6 loop mutation is easily under¬ 
stood from the structure of the BCoV HE-receptor complex, as 
it involved loss of Y 184 , an essential residue for ligand binding. 

Crystal Structure of an 0C43A-type HE 

To assess the consequences of the combined mutations in the 
type A HE lectin domain, we determined the crystal structure of 


OC43 NL/A/2005 HE at 2.45A resolution 
(for crystallographic details, see Fable 
SI; PDB: 5N11). In overall fold, OC43 
HE is highly similar to other coronavirus 
HEs ( : igure 4). Differences in the 
esterase domain are limited to substitu¬ 
tion of five surface residues, none likely 
to impair esterase function or activity. 
OC43 HE is also identical to BCoV HE 
in most of its lectin domain, particularly 
in the beta-sandwich core, but the loops 
shaping the receptor-binding site have changed dramatically 
( igure 4D; r igure S3D). The pi 1-pi 2 and p5-p6 loops are 
disordered, suggestive of extensive flexibility, in the case of 
the latter loop in accordance with loss of E 178 -mediated inter¬ 
loop hydrogen bonding. Those (parts of the) loops that can be 
visualized are reoriented and displaced with respect to their 
original position in BCoV HE, utterly destroying the CBS. For 
example the p7/p10 p-hairpin, comprising F 211 , L 212 , S 213 , and 
N 214 that are key to protein-ligand interaction in BCoV HE 

o 

( !eng et al., 2008), has shifted by approximately 7 A ( igure 4D). 
The PI pocket, arguably the most critical element of the CBS, 
no longer exists as (1) F 211 is no longer at its original position 
and its side chain is rotated by 90°, (2) Y 184 was lost as a result 
of the frameshift mutation in the p5-p6 loop, and (3) the side 
chains of L 266 and L 267 are reoriented apparently due to the 
loss of the MBS ( igures 4E and 4F). Finally, attached to N 114 
is an A/-acetylglucosamine, providing formal evidence that the 
newly introduced glycosylation site in OC43 HEs is functional. 
Although the remaining residues of the sugar chain could not 
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be visualized, the glycan would stretch across the remains of 
the CBS (Figure 4F). 

Loss of HE Receptor Binding in HCoV-HKUl: A Case of 
Convergent Evolution? 

OC43 and HKU1 occupy the same niche, prompting the ques¬ 
tion whether they have been subject to similar selective con¬ 
straints. HCoV-HKUl is also a lineage A betacoronavirus but 
distantly enough related to the (31 CoVs to be assigned to a sepa¬ 
rate species (de Groot et al., 2011) ( : igure 1A). Like OC43, 
HCoV-HKUl binds to 9-O-Ac-Sia receptor determinants via its 
spike (Huang et al., 2015), and it possesses an HE (Woo et al., 
2005b). The HE of the prototypical HCoV-HKUl 2005 Hong 
Kong genotype A strain (Woo et al., 2006) tested negative for 
9-O-Ac-Sia binding by sp-LBA, standard HAA, and nanobead 
HAA ( igures 1 B-1D). Comparative sequence analysis, including 
genotypes A and B, of HCoV-HKUl (Woo et al., 2006) shows that 
the HKU1 HE has undergone massive deletions, as a result of 
which most of the lectin domain was lost ( -igure S4). 

Loss of HE Lectin Function Reduces Receptor- 
Destroying Enzyme Activity toward Multivalent 
Substrates 

In influenza C and D viruses, the hemagglutinin-esterase fusion 
protein is uniquely responsible for viral attachment and entry, 
and hence HE function is receptor binding first and foremost. 
The same holds for the HEs of MHV-type murine betacoronavi- 
ruses (MuCoV) that mediate virion binding to O-Ac-sialoglycan- 
based attachment factors (Langereis et al., 2010), while S binds 
to the proteinaceous entry receptor CEACAM-1 (Peng et al., 
2011). At the other end of the spectrum, in OC43 and HKU1, 
HE’s sole remaining function is receptor destruction, with virion 
attachment to 9-O-Ac-sialoglycans exclusively assigned to S. 
For the animal (31 CoVs that retain a functional HE lectin domain, 
however, the role of HE is less obvious. On the one hand, HE 
might participate in attachment to increase virion binding avidity. 
On the other, its lectin domain may serve primarily to promote 
catalytic activity toward high-multivalency substrates by 
bringing the esterase in close proximity to and prolonged asso¬ 
ciation with clustered glycotopes, in analogy with the carbohy¬ 
drate-binding modules of cellular glycoside hydrolases (Boras- 
ton et al., 2004). In the latter case, the HE lectin CBS in the 
human lineage A (3CoVs may have been lost to downregulate re¬ 
ceptor-destroying activity. To study this, we compared soluble 
recombinant HEs of BCoV, OC43, and HKU1 for their relative 
enzymatic activity toward the monovalent substrate p-nitrophe- 
nol acetate (pNPA) and the multivalent substrate BSM. The latter 
glycoconjugate carries hundreds of O-linked glycans clustered 
in such densities to give the protein a filamentous bottlebrush 
appearance typical for mucins (lappone et al., 2015). Whereas 
all three HEs showed comparable specific activities when as¬ 
sayed with pNPA ( : igure 5A), those of OC43 and HKU1 were 
more than 250-fold less active than BCoV HE in receptor- 
destruction assays with BSM ( Igure 5B). Inactivation of the 
BCoV HE lectin CBS through a F 211 A substitution reduced its 
esterase activity toward BSM to that of wild-type OC43 HE. 
Conversely, restoration of the OC43 HE lectin domain increased 
esterase activity to that of BCoV HE (see “TREF,” Figure 5B). 
Thus, BCoV and HCoV HEs indeed differ in their reactivity toward 


clustered receptors in consequence of the respective conserva¬ 
tion or loss of lectin function. 

Consequences of Loss of Lectin Function for Virion- 
Mediated Receptor Destruction 

OC43 HE efficiently depletes endogenous pools of O-Ac-Sias 
within infected cells ( igure S5) as is critical for the release of viral 
progeny (Desforges et al., 2013). The loss of HE lectin function 
therefore is more likely to affect virus biology at another stage, 
namely that of virion (pre)attachment. Of note, the observations 
for recombinant soluble HEs cannot be extrapolated one-to- 
one to virion-associated HEs, because in virus particles, HE 
esterase activity will be affected by S-mediated attachment to 
9-O-Ac-Sias as well as by the close packing of HE molecules 
in the viral envelope. Serial dilutions of intact BCoV and OC43 vi¬ 
rions, adjusted for pNPA esterase activity (Table S2), were there¬ 
fore compared in solid-phase receptor-destruction assays with 
BSM ( : igures 5C and 5D). Remarkably, the difference between 
the two viruses as measured by the onset of detectable receptor 
destruction was smaller than for their respective soluble HEs. 

The activity of virion-associated BCoV HE relative to that of 
soluble BCoV HE was reduced 8- to 16-fold. Conversely, the 
relative activity of virion-associated OC43 was increased at least 
in that destruction of receptors already became apparent at 
lower enzyme concentrations. ( : igure 5C). To study whether 
these phenomena relate to the S protein and to S-mediated 
virion attachment to 9-O-Ac-Sias, we constructed recombinant 
murine coronaviruses (rMuCoVs) that express MuCoV S, which 
does not bind 9-O-Ac-Sias, in combination with either BCoV or 
OC43 HE. rMuCoV virions studded with OC43 HE lacked detect¬ 
able esterase activity toward BSM. In contrast, the esterase ac¬ 
tivity of rMuCoVs expressing BCoV HE was 6- to 8-fold higher 
than that of wild-type BCoV virions ( igure 5C). Apparently, 
|31CoV S proteins, by binding to 9-O-Ac-Sias negatively or posi¬ 
tively, affect HE esterase activity in BCoV and OC43 virions, 
respectively, depending on whether or not the HE lectin domain 
is functional. Most strikingly, BCoV and OC43 virions differed in 
the extent to which receptors were depleted. The soluble HEs of 
OC43 and BCoV, despite their difference in enzyme activity, de¬ 
stroyed all O-Ac-Sias as detectable by sp-LBA. BCoV virions de¬ 
stroyed 95% of the receptors after 24 hr incubation, which 
increased to 98% after 48 hr ( igure 5D). For OC43 virions, the 
effect was more pronounced such that after 24 hr only 85% of 
the receptors were destroyed. The remaining receptors were 
resistant even to continued incubation ( : igure 5D). The data sug¬ 
gest that, in the context of the virion, the absence of a functional 
HE lectin domain results not only in reduced enzyme activity to¬ 
ward clustered receptors, but also prohibits cleavage of partic¬ 
ular receptor populations. As these “protected” receptor popu¬ 
lations are readily depleted by soluble OC43 HE ( : igures 5B-5D), 
their resistance to cleavage by OC43 virions suggests that they 
are inaccessible to HEs when embedded in the viral envelope 
and in the presence of S. 

The data prompt a model in which the susceptibility of popu¬ 
lations of clustered receptors to cleavage depends on (1) their 
accessibility, as determined by the distance by which they 
extend from a fixed surface—the bottom of the ELISA well in 
our artificial system and, for example, the host cell membrane 
in vivo—and (2) the considerable difference in height between 
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Figure 5. Loss of HE Lectin-Mediated 
Receptor Binding Alters Sialate-9-O-Acety- 
lesterase Receptor, Destroying Activity 
toward Multivalent Substrates 

(A) Esterase activity of soluble recombinant 
HE + -Fc fusion proteins (sHE + ) toward monovalent 
substrate pNPA. Enzymatic activity is shown as 
percentage of BCoV HE wild-type activity. 
WT, sHE + with wild-type HE ectodomains of 
BCoV-Mebus, OC43 USA/1967, or HCoV-HKUl 
as indicated; F 211 A, BCoV-Mebus sHE + derivative 
with the lectin CBS inactivated through a 
F 211 A substitution ( !eng et al., 2008); TREF, 
OC43 USA/1967 HE with the lectin CBS repaired 
(as in Figure 2C). Data are represented as 
mean ± SEM. 

(B) Loss of sHE + -mediated receptor binding re¬ 
sults in reduced esterase activity toward high- 
multivalency substrates as determined by on-the- 
plate O-Ac-Sia depletion assay. BSM, coated on 
maxisorp plates, was incubated for 1 hr with 
2-fold serial dilutions of sHE + . Receptor destruc¬ 
tion was assessed by sp-LBA with BCoV HE°-Fc 
at fixed concentration (2 ng/^iL). 

(C) Whole-virion receptor-destruction assays 
with 2-fold serial dilutions of purified viruses 
(starting at 15 mU). BCoV and OC43 sHE + s were 
included for comparison. Receptor destruction 
was measured after 1 hr or 4 hr incubation. Black 
dotted lines indicate 100% receptor destruction 
as determined with excess amounts of BCoV or 
OC43 sHE + . 

(D) Whole-virus-mediated receptor destruction 
over time with fixed amounts of BCoV-Mebus and 
OC43 USA/1967 virions (15 mU). To exclude 
“exhaustion” (inactivation of OC43 esterase 
and/or changes in the physicochemical properties 
of the virions over time), we removed BCoV and 
OC43 at t = 24 and replaced them with equal 
amounts of freshly thawed aliquots of the virus 
stocks (black arrow). 

(E) Scaled side-by-side representations of the 
structures of coronavirus S (Walls et al., 201 6) and 
HE proteins. Indicated are the estimated heights 
of S and HE and the approximate distance sepa¬ 
rating S receptor-binding sites and HE catalytic 
sites. (See also Figure S4 and S5 and Table S2.) 


S and HE (20 and 8 nm, respectively), with the CBSs of S and 
those of HE esterase domains separated by a distance of up 
to 15 nm ( igure 5E). S-mediated virion attachment to a 
9-O-Ac-sialylated surface would bring the HE catalytic site in 
proximity of some clustered receptors (which would explain 
the stimulating effect of S on HE esterase activity in the case 
of OC43 virions, but not observed for rMuCoV-OC43 HE + ) but 
at the same time would keep the virus-associated enzyme at 
“arm’s length” from other receptor populations, readily acces¬ 
sible to soluble HEs (which would explain the inhibiting effect 
of S on esterase activity in the case of BCoV virions and the 
apparent lack thereof in the case of rMuCoV-BCoV HE + ) ( igures 
5C and 5D). Thus, a distinction arises between clustered recep¬ 
tor populations that come within immediate reach of the HE 
esterase and that therefore would be cleaved by virion-associ¬ 
ated HE as efficiently as by soluble HE and those that are kept 
out of reach of virion-associated HE as a result of S-mediated 


receptor-binding and that are therefore cleaved at strongly 
reduced rates or even rendered non-cleavable unless HE is pro¬ 
vided with a functional lectin domain ( : igure 6 ). By virions 
grasping on to multivalent sialoglycoconjugates via the HE re¬ 
ceptor-binding sites and “drawing them in” closer to the surface 
of the envelope, clustered receptors would yet become available 
to the HE esterase domain ( igure 6 ). Such a mechanism may be 
promoted by cooperativity among adjacent HE molecules and, 
membrane fluidity permitting, possibly even lead to redistribu¬ 
tion of surface projections in the viral envelope, resulting in 
displacement of spikes, local recruitment of HEs, and formation 
of a “receptor-destruction patch.” Saliently, some receptor pop¬ 
ulations may be completely resistant to destruction even by 
virion-associated HEs with a functional lectin domain, such as 
the 2% of detectable O-Ac-Sias remaining after 48 hr of incuba¬ 
tion with BCoV ( igure 5D). It is tempting to speculate that under 
natural conditions, the distinction between decoy and entry 
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Figure 6. Hypothetical Model for the Inter¬ 
action of HCoV and BCoV Virions with Multi¬ 
valent Glycoconjugates 

Schematically depicted are portions of HCoV 
(OC43 or HKU1) and BCoV virions, with large 
spikes comprised of S (in gray) and smaller pro¬ 
trusions comprised of HE (yellow) extending from 
the viral membrane (blue). Functional HE lectin 
CBSs are indicated by black holes (only one 
shown per HE dimer for reasons of simplicity). Also 
shown schematically are membrane-anchored 
(bottom) and non-anchored (top) mucin-type 
glycoconjugates of bottle-brush filamentous 
appearance with clustered receptors (9-O-Ac- 
Sias, red dots) arranged in linear arrays and with 
absence of red dots indicating receptor-destruc¬ 
tion. The model, based on the size difference be¬ 
tween S and HE, visualizes how loss of HE lectin 
function might alter virion-associated receptor- 
destroying activity, reducing the specific activity of 
virions as well as the rate and selectivity of re¬ 
ceptor destruction. In virions with lectin-deficient 
HEs, clustered substrates will be largely kept at a 

distance from the HE esterase catalytic pocket as a result of S-glycoconjugate interaction. In contrast, HEs with intact lectin CBS may draw in portions of the 
glycoconjugates (or will draw the virion-associated HEs toward them), aided by cooperativity of binding between adjacent HEs within the viral envelope. Thus, 
clustered glycotopes become fixed within reach of the esterase catalytic sites and receptor destruction is promoted. 


receptors is made by default on basis of their accessibility to 
cleavage by virion-associated HE. 

CONCLUSION 

We showed that OC43—after its zoonotic introduction —has 
been under incessant selective pressure to adapt to the sialogly- 
come of the human host. HE-associated receptor binding was 
selected against, largely lost early on, and ultimately lost alto¬ 
gether through an accumulation of mutations in the lectin domain 
over a period of decades. In result, the balance between virion 
attachment and catalytically driven virion elution was reset, 
apparently to meet the specific requirements for optimal replica¬ 
tion in human airways and to allow the virus in this particular 
niche to distinguish between decoys marked for destruction 
and functional receptors that should be preserved for cell entry. 
This view is reinforced by our observation that HCoV-HKUl, 
another respiratory lineage A pCoV of humans, followed a 
convergent evolutionary path and also lost HE-mediated recep¬ 
tor binding. We offer that differences in the sialoglycomes of bo- 
vids and humans, such as in the local expression, structure, and 
density of 9-O-acetylated sialoglycoconjugates, may pose in¬ 
compatibilities during cross-species transmission of OC43 and 
BCoV and that the respective loss or preservation of HE-medi- 
ated receptor binding contributes to the host tropism of these 
viruses. This said, comparative sialoglycomics is still in its in¬ 
fancy, humans and bovids have not been studied for differences 
in O-Ac-sialoglycan expression in any detail, and the precise sia- 
loglycan-based constraints that selected for the particular traits 
of the animal and human pCoVs therefore remain to be identified. 
At any rate, our findings do reveal an as yet unappreciated 
aspect of lineage A pCoV adaptation to humans. Of broader 
relevance, they provide a general paradigm for dynamic, catal¬ 
ysis-driven virus-sialoglycan interactions in relation to host 
selectivity. This notion is supported by observations for influenza 


A viruses, where neuraminidase (NA) activity toward multivalent 
substrates varies with differences in the length of the stalk 
domain and positively correlates with NA size (Castrucci and 
Kawaoka, 1993; Els et al., 1985), where NA catalytic activity is 
modulated by the absence or presence of a second Sia-binding 
site (Uhlendorff et al., 2009; Varghese et al., 1997) and where 
changes in NA stalk length and NA lectin affinity have been impli¬ 
cated in host specificity (Blumenkrantz et al., 2013; Castrucci 
and Kawaoka, 1993; Lai et al., 2012; Li et al., 2014; Matrosovich 
et al., 1999; Munier et al., 2010; Uhlendorff et al., 2009; Varghese 
et al., 199' ). These observations can be readily interpreted in the 
context of our data and model proposed. 

EXPERIMENTAL PROCEDURES 

Materials and experimental procedures are detailed in the Supplemental 
Experimental Procedures. 

Expression and Purification of HEs 

Human codon-optimized sequences of the HEs of BCoV-Mebus, ECoV-NC99, 
PHEV-VW572, RbCoV-HKU14-1, CRCoV-240/05, HCoV-HKUl, and OC43 
USA/1967 (ATCC VR-759) were cloned in expression plasmid pCD5-T-Fc 
(7eng et al., 2008). The resulting constructs encode chimeric proteins 
comprised of the HE ectodomain fused to the human IgGI Fc domain, with 
the domains separated by a thrombin cleavage site. The fusion proteins 
were either expressed in an enzymatically active form (HE + ) or rendered inac¬ 
tive through a catalytic Ser-to-Ala substitution (HE 0 ). Plasmids coding for 
OC43 NL/A/2005 and OC43 NL/B/2005 HE were constructed through site- 
directed mutagenesis of the OC43 USA/1967 HE expression vector using 
the Q5 kit (New England Biolabs). For lectin affinity and esterase activity as¬ 
says, HE-Fcs were produced by transient expression in HEK293T cells and pu¬ 
rified from cell culture supernatants by protein A affinity chromatography and 
low pH elution as described (.eng et al., 2008). For crystallization, OC43 
NL/A/2005 HE + -Fc was expressed in HEK293 GnTI(-) cells and purified by 
protein A affinity chromatography, followed by on-the-beads thrombin cleav¬ 
age as described (Zeng et al., 2008). Beads were pelleted, and the HE ectodo¬ 
main in the supernatant was concentrated to ^15 mg/mL. HE genes from 
OC43 field strains were RT-PCR amplified with OC43 viral RNA, directly 
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isolated from nasopharyngeal aspirates (van de Pol et al., 2006) as a template, 
and their sequences were used to construct codon-optimized versions. 

Purification of Virions for Whole-Virus Receptor-Destruction Assays 

Cell monolayers were infected at MOI 0.01. Virions were purified by ultracen¬ 
trifugation through 20% (w/v) sucrose cushions and resuspended in PBS. 
Virus preparations were analyzed for particle content by qPCR, plaque 
assay, quantitative latex bead ratio EM, and pNPA esterase activity assay 

(Table S2). 

Hemagglutination Assay 

HAA was performed with rat erythrocytes (Rattus norvegicus strain Wistar; 
50% suspension in PBS) and, in standard tests, with 2-fold serial dilutions of 
HE°-Fc proteins (starting at 25 ng/well) as described ( leng et al., 2008). For 
high-sensitivity HAA, based on multivalency-driven high-avidity binding, 5 jig 
HE°-Fc was complexed to 2 x 10 9 protein A-coated 100 nm nanobeads 
(Chemicell) in 100 jil_ PBS for 45 min at 4°C prior to 2-fold serial dilution. Hem¬ 
agglutination was for 2 hr at 4°C unless stated otherwise. 

Solid-Phase Lectin-Binding Assay and On-the-Plate O-Ac-Sia 
Depletion Assay 

sp-LBAs were performed with bovine submaxillary mucin (BSM) and 2-fold se¬ 
rial dilutions of HE°-Fc proteins as described (Langereis et al., 2012, 2015). 

Receptor-destroying esterase activities of soluble recombinant HEs (sHE+) 
and whole viruses were measured by on-the-plate O-Ac-Sia depletion assay, 
performed essentially as described (Langereis et al., 2015). BSM, coated on 
Maxisorp flat-bottom 96-well plates, was (mock-)treated with 2-fold serial di¬ 
lutions of sHE + in PBS (100 jiL/well) starting at 20 ng/jiL for 1 hr, or, for whole- 
virion receptor-destruction assays, with virions or sHE + starting at 15 mU 
esterase activity for up to 48 hr at 37°C. Depletion of O-Ac-Sia receptors 
was detected by sp-LBA with BCoV-LUN HE 0 (2 ng/|iL; 100 (.iL/well). Esterase 
activity was plotted as the inverse of lectin-binding measurements, expressed 
in percentages. 
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